Approximate Inverses and Incomplete Decompositions
نویسنده
چکیده
In this work we consider the approximate inverse preconditioners AISM based on the Sherman-Morrison formula [1]. We will show some of theoretical properties of the factors of AISM in function of the value of parameter s considered in the decomposition. In particular, we will show that one of the factors hides a scaled factor of the standard incomplete LDU decomposition of the system matrix A. Thus the AISM decomposition can be considered as a way to get an incomplete decomposition via a factorizad approximate inverse. We will show that the LDU decomposition obtained in this way can be useful in preconditioned iterative methods. Numerical experiments illustrate the efficiency of this incomplete decomposition where some dropping strategies are used where the entries and the norms of the inverse factors are considered. Acknowledgement: This research was fully supported by Spanish grantDGI grant MTM2004-02998 and by the project 1ET400300415 within the Na-tional Program of Research “Information Society” References[1] R. Bru, J. Cerdn, J. Marn and J. Mas, Preconditioning nonsymmetric linearsystems with the Sherman-Morrison formula, SIAM J. Sci. Comput., 25, 2 (2003),pp. 701–715.
منابع مشابه
Sparse Approximate Inverses for Preconditioning of Linear Equations
1. Sparse Approximate Inverses and Linear Equations We consider the problem of solving a system of linear equations Ax = b in a parallel environment. Here, the n n-matrix A is large, sparse, unstructured, nonsymmetric, and ill-conditioned. The solution method should be robust, easy to parallelize, and applicable as a black box solver. Direct solution methods like the Gaussian Elimination are no...
متن کاملApproximate Inverse Preconditioners via Sparse-Sparse Iterations
The standard incomplete LU (ILU) preconditioners often fail for general sparse in-deenite matrices because they give rise tòunstable' factors L and U. In such cases, it may be attractive to approximate the inverse of the matrix directly. This paper focuses on approximate inverse preconditioners based on minimizing kI ? AMk F , where AM is the preconditioned matrix. An iterative descent-type met...
متن کاملOn the Relations between ILUs and Factored Approximate Inverses
This paper discusses some relationships between Incomplete LU (ILU) factoriza-tion techniques and factored sparse approximate inverse (AINV) techniques. While ILU factorizations compute approximate LU factors of the coeecient matrix A, AINV techniques aim at building triangular matrices Z and W such that W > AZ is approximately diagonal. The paper shows that certain forms of approximate inverse...
متن کاملOn the numerical solution of convection-dominated problems using hierarchical matrices
The aim of this article is to shows that hierarchical matrices (H-matrices) provide a means to efficiently precondition linear systems arising from the streamline diffusion finite-element method applied to convection-dominated problems. Approximate inverses and approximate LU decompositions can be computed with logarithmic-linear complexity in the standard Hmatrix format. Neither the complexity...
متن کاملA Comparative Study of Sparse Approximate Inverse Preconditioners
A number of recently proposed preconditioning techniques based on sparse approximate inverses are considered. A description of the preconditioners is given, and the results of an experimental comparison performed on a Cray C98 vector processor using sparse matrices from a variety of applications are presented. A comparison with more standard preconditioning techniques, such as incomplete factor...
متن کامل